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In this paper we study the double J/iji production in central diffractive processes considering the 
Resolved Pomeron model. Based on the nonrelativistic QCD (NRQCD) factorization formalism for 
the quarkonium production mechanism we estimate the rapidity and transverse momentum depen¬ 
dence of the cross section for the double J/ij) production in diffractive processes at LHC energies. 
The contributions of the color-singlet and color-octet channels are estimated and predictions for 
the total cross sections in the kinematical regions of the LHC experiments are also presented. Our 
results demonstrate that the contribution of central diffractive processes is not negligible and that 
its study can be useful to test the Resolved Pomeron model. 
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The study of the production of heavy quarkonium states provides a unique laboratory in which one can explore 
the interplay between perturbative and nonperturbative effects in QCD (For a review see, e.g., Ref. i). In the 
last years, several theoretical approaches have been proposed for the calculation of these states, as for instance, the 
Color Singlet model (CSM), the color evaporation model, the Non Relativistic QCD (NRQCD) approach, the frag¬ 
mentation approach, and the fey-factorization approach. Albeit considerable efforts both in theory and experiments, 
the quarkonium production mechanism is still not fully understood. A particular example is the double quarkonium 
production in inclusive processes, i. e. processes where the two incident hadrons dissociate in the interaction. In 
the last years, the measurements reported by the LHCb Q, CMS Q and DO Q collaborations at the LHC and the 
Tevatron have posed significant challenges to our understanding of the quarkonium production. Currently, there are 
in the literature several predictions for the total cross section and differential distributions [Ml, at leading-order 
and next-to-leading order of the perturbative expansion, some of them using the Color Singlet model and others using 
the NRQCD formalism. The results of these studies indicate that the discrimination between the different approaches 
in inclusive processes will be a hard task, since double parton scattering (DPS) processes are expected to contribute 
at high energies and the contribution of this new mechanism is still an open question [Tsl - I^ . In contrast, the DPS 
contribution for diffractive processes, where the incident hadrons remain intact, is expected to be negligible [^ . This 
motivates the study of the double quarkonium production in diffractive interactions in order to test the production 
mechanism. 

In recent years the diffractive processes have attracted much attention as a way of amplif ying the physics programme 
at hadronic colliders, including searching for New Physics (For a review see, e.g. Ref. [^). The investigation of 
these reactions at high energies gives important information about the structure of hadrons and their interaction 
mechanisms. The diffractive physics has been tested in hadron-hadron collisions considering different theoretical 
approaches and distinct final states like dijets, electroweak vector bosons, dileptons, heavy quarks, quarkonium -|- 
photon and jet -|- photon (See, e.g.. Refs . One of these approaches is the Resolved Pomeron Model, 

proposed by Ingelman and Schlein in Ref. [39j, which assumes the validity of the diffractive factorization formalism 
and that the Pomeron has a partonic structure. The basic idea is that the hard scattering resolves the quark and 
gluon content in the Pomeron which can be obtained by analysing the experimental data from diffractive deep 
inelastic scattering (DDIS) at HERA, providing us with the diffractive distributions of singlet quarks and gluons in 
the Pomeron (^ . However, other approaches based on very distinct assumptions are also able to describe the current 
scarce experimental data. Consequently, the present scenario for diffractive processes is unclear, motivating the study 
of alternative processes which could allow to constrain the correct description of the Pomeron. 

In this paper we propose the study of the double J/T* production as a complementary test of diffractive processes 
and the pomeron structure. In particular, we present a detailed analysis of the rapidity and transverse momentum 
dependence of the cross section for the double J/Al production in diffractive processes considering the nonrelativistic 
QCD (NRQCD) factorization formalism [4l| for the quarkonium production mechanism. Moreover, a comparison 
with the inclusive production is presented. Our analysis is strongly motivated by the recent LHCb data for the 
diffractive production of this final state and by the recent theoretical study performed in Ref. [2l| , which have studied 
the double J/Al production in exclusive processes. Our goal is to present a complementary analysis using a distinct 
model for the treatment of the diffractive interactions. 

In the following we apply the Resolved Pomeron Model for the central diffractive double J/ij} production. 
This model assumes that the Pomeron has a well defined partonic structure and that the hard process takes place 
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FIG. 1: Typical diagrams for the double J/ijj production in (a) central diffractive and (b) inclusive processes. 


in Pomeron - Pomeron processes. The contributing diagrams are the same as in the inclusive case, and the J/ip 
production in the central diffractive processes is described by diagrams like those ilustrated in Fig. [T](a), with IP IP 
interactions. Differently from the inclusive processes illustrated in Fig. [T](b), the diffractive processes are characterized 
by the presence of one or two intact very forward hadrons and empty regions in pseudo-rapidity, called rapidity gaps, 
in the final state. In this paper we restrict our study to the central diffractive processes, which are characterized by 
two intact hadrons, two rapidity gaps and soft particles accompanying the double J/'i' production [See Fig. [T](a)]. 
The presence of these soft particles associated to the Pomeron remnants is a characteristic of the Resolved Pomeron 
model. In contrast, in the case of the central exclusive production of two J/T' discussed in Ref. [2l|, nothing else is 
produced except the leading hadrons and the central object. 

In the NRQCD formalism (dTj the cross section for the diffractive production of a heavy quarkonium pair can be 
factorized as follows 

da{pp -)-p'^Xi+Hi+H 2 + X 2 'Sip) = {S'^) ^ f^/p ® f^/p 0 da[ab -)■ QQ[ni] + QQ[n 2 ] + X] ■ ,(1) 

a,6,ni ,n2 


where ® represents the presence of a rapidity gap in the final state, Xi the remnants of the Pomeron, (5^) is 
the gap survival probability (see below), are the diffractive parton distributions and the coefficients da\ab —>■ 
QQ[ni] + QQ[n 2 ] + X] are perturbatively calculable short distance cross sections for the production of the two heavy 
quark pairs in an intermediate Fock state rii = _ i^g)^ which does not have to be color neutral. The 

{O^') are nonperturbative long distance matrix elements (LDME), which describe the transition of the intermediate 
QQ in the physical state H via soft gluon radiation. Currently, these elements have to be extracted from global fits 
to quarkonium data as performed, for instance, in Ref. [i^. In the Color Singlet Model [ij, only the states with the 
same quantum numbers as the resulting quarkonium contribute to the formation of a bound QQ state. In contrast, in 
NRQCD, also color octet QQ states contribute to the quarkonium production cross section via soft gluon radiation. 
The Color Singlet Model can be obtained from NRQCD factorization by retaining, for a given process, only the 
contribution that is associated with the color-singlet LDME of the lowest non-trivial order in v, which is the typical 
velocity of the heavy quark or antiquark in the quarkonium rest frame. At high energies the double quarkonium 
production is dominated by gluon - gluon interactions. Consequently, the differential cross section for the double J/^/> 
production in central diffractive (CD) processes can be written as 


da^° 

dydp^ 



dXaQ {Xa-^P ')p (^ 6 ,/^ ) 


D/ 


XaXb 


2xa — XTe.y 


^ ^^[gg^2cU^SQ\-{Oi/^eS,))\ 


( 2 ) 


where Xamin = 2 -sfe-y ’ = 2 xa-Lly ’ ^ +Pt- Here M is the J/ip mass, px its 

transverse momentum and y its rapidity. The J/')/’ transverse mass is taken as the hard scale of the problem, with 
Pr = PF = rriT- Here, g^{xi,p^) are the diffractive gluon distribution functions from the two colliding protons. 
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In the present work, the diffractive gluon distributions in the proton are taken from the Resolved Pomeron Model 
[ 3 ^, where they are defined as a convolution of the Pomeron flux emitted by the proton, fpixp), and the gluon 
distribution in the Pomeron, gp{l3, where /3 is the momentum fraction carried by the partons inside the Pomeron. 
The Pomeron flux is given by fp{xp) = dtfp/p{xp,t), where fp/p{xp,t) = Ap ■ and tmin, tmax are 

■min Xjp 

kinematic boundaries. The Pomeron flux factor is motivated by Regge theory, where the Pomeron trajectory assumed 
to be linear, ap{t) = ap{0) + a'pt^ and the parameters Bp, a'p and their uncertainties are obtained from fits to HI 
data (40j| . The diffractive gluon distribution is then given by 


g^{x,g^) = [ dxpd(35{x- xpP)fp{xp)gp{l3,g^) = [ ^^^fp{xp)gp( — ,nA . (3) 

J Jx Xp \xp J 

In our analysis we use the diffractive gluon distribution obtained by the HI Collaboration at DESY-HERA [i^. It 
is important to emphasize that the cross section for the diffractive production of two J/'k is strongly sensitive to the 
Pomeron structure due to the quadratic dependence on the diffractive gluon distribution. Finally, ^ in Eq. © are 
the hard scattering differential cross sections, which we assume to be given by leading order (LO) af expressions. 
As demonstrated in Ref. the corresponding Feynman diagrams can be classified into two groups: (a) diagrams 
associated to the nonfragmentation contribution, composed by 31 Feynman diagrams, with the leading contribution 
being the color singlet (cc)i(^S'i) + (cc)i(^«S'i) channel, and (b) diagrams associated to the gluon fragmentation 
contribution, composed by 72 Feynman diagrams, with the main contribution associated to the color octet ( 00 ) 8 (^ 51 ) + 
(cc)8(^<S'i) channel. For the gluon-initiated color singlet contributions, one has 0] 


^[gg^2cc,es^)]■{o(/^es,r = 

at 


167rag|i?(0)|^ 

81M2s8(M2 - t)4(M2 - m)4 


jki 


( 4 ) 


where, as in the Ref. Q, |i?(0)|2 = 0.8 GeV^ is the squared radial function at the origin, s, t and u are the usual 
Mandelstam variables, M = 2toc, and rric = 1.5 GeV. The detailed expressions for the ajki coefficients in Eq. ([4]) 
are given in Ref. Q- On the other hand, for the color octet contributions, the differential cross section for the gluon 
initiated partonic snbprocesses was calculated in Ref. m and can be written as 


da 

di 


[gg -)> 2cc8('^S'i)] = 




972M658(M2 - i)4(M2 - w)4 


14 

E 

j=0 




( 5 ) 


where the aj coefficients can be found in Ref. M- In the considered case, the only relevant NRQCD matrix element 
is (Og^’^(^S'i)) = 3.9 X 10“^GeH^, taken from [i^. This value has been updated in a recent fit to world data 
which gives = 1.68 x 10“^GeF^. Using this new value, our results for the COM contributions would 

decrease by a factor 2.3, which is inside the uncertainties of our results. Moreover, it is important to emphasize 
that, as in Ref. 0, we have disregarded the contributions from the So)) transition. This approximation is 

reasonable at high transverse momentum pT- In contrast, for small values of pT, these contributions are not negligible, 
dominating the COM contributions, as demonstrated in Refs. Consequently, our COM predictions, denoted 

COM* in what follows, shonld be considered an incomplete evaluation of the COM contributions at \ow-pT- However, 
it is important to emphasize that the inclusion of the (Og^’^(^S'o)) transition in our calculations would not modify our 
main conclusions, since the CSM contribution largely dominates the cross section at small values of the transverse 
momentum (see below). 

In order to obtain realistic predictions for the double J/T production in central diffractive processes, it is crucial to 
use an adequate value for the gap survival probability, (5^). This factor is the probability that secondaries, which are 
produced by soft rescatterings do not populate the rapidity gaps, and depends on the particles involved in the process 
and in the center-of-mass energy. In what follows, following Ref. [i^, we will assume that (5^) = 2% for proton - 
proton collisions at LHC energies. However, this snbject deserves a more detailed analysis, since the magnitude of 
( 52 ) is still a theme of intense debate in the literature (See, e.g.. Ref. [47j|l. 

In Fig. [2] we present our results for the transverse momentum distribution for the double J/'k production at 
midrapidities (|?/| < 2.5) in central diffractive processes for LHC energies. For the sake of comparison, we also show 
the results for the inclusive case, obtained using the CTEQ6L parametrization for the gluon distribution in the 
proton, which agree with the predictions presented in Ref. 0- The central diffractive predictions at small px are a 
factor Ri 10^ smaller than in the inclusive one. We obtain that the px distributions in the low px region are dominated 
by the color singlet contributions. Moreover, as in the inclusive case 0, the distribution vanishes at px = 0 and 
increases rapidly until it reaches a maximum at px ~ 1-5 GeV. Then it decreases monotonically as px increases. We 
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Central diffractive double J/\|/ @ LHC Central diffractive double J/\|r @ LHC 

8 TeV, lyl<2.5 14 TeV, lyl<2.5 



FIG. 2: (Color online) Transverse momentum distributions for double J/ip production in central diffractive processes at -y/s = 8 
TeV (left panel) and 14 TeV (right panel). The prediction for the inclusive production is presented for comparison. 
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FIG. 3: (Color online) Rapidity distributions for double J/'ip production in inclusive and central diffractive processes in pp 
collisions at ^/s = 8 and 14 TeV. 


obtain that at leading order the color singlet contributions are dominant except at large pp- The cross-over, beyond 
that the color octet contributions start to be dominant, occurs at pp ~ 15 GeV. Furthermore, in the dominant \ow-pp 
peak, the color octet contribution is four orders of magnitude less important than the color singlet one. It is important 
to emphasize that the recent studies of the next-to-leading order corrections for the inclusive double J/T' production 
performed in Ref. [lB | indicate the the color singlet contributions also are dominant at large-pr- A similar behaviour 
also is expected in the case of the double J/T' production in central diffractive processes if the NLO corrections are 
taken into account. 

In Fig. 13] we present our results for the rapidity distribution for the double J/4> production in central diffractive 
processes at \/s = 8 and 14 TeV. The inclusive predictions are also presented for comparison. As expected our 
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CM energy 

Kinematical ranges 

Inclusive 

Central Diffractive 

CEP 

8 TeV 

\v\ < 2.5 

27203 pb 

9.51 pb 

10 pb 

8 TeV 

|y| < 2.5 and pT > 6.0 GeV 

29 pb 

4.5 X 10“"’ pb 

- 

8 TeV 

2 < y < 4.5 

9709 pb 

2.16 pb 

2.5 pb 

14 TeV 

\V\ < 2.5 

39690 pb 

16.02 pb 

17 pb 

14 TeV 

|y| < 2.5 and pT > 6.0 GeV 

46 pb 

9.2 X lO^"* pb 

- 

14 TeV 

2 < y < 4.5 

15220 pb 

4.43 pb 

4.7 pb 


TABLE I: Total cross sections for double J/xb production for different energies and distinct rapidity and transverse momentum 
cuts. Also shown the predictions from Ref. for the central exclusive production (CEP). 


predictions increase with the energy. We have verified that the color singlet contributions dominate all regions for 
the pT-integrated spectra, with the color octet one being negligible in all rapidity regions. Moreover, we have a 
reduction of four orders of magnitude when going from the inclusive to the central diffractive case. Our results for 
the rapidity distribution allows us to obtain the predictions for the cross section in distinct rapidity ranges covered 
by the different LHC experiments. In Table U we present our predictions for the total cross section considering the 
production at midrapidities (|?/| < 2.5) , which can be analyzed by CMS, ATLAS and ALICE Collaborations, and for 
the kinematical range probed by the LHCb Collaboration {2 < y < 4.5). For comparison we present the predictions 
for the inclusive production and for the double J/T' production in central exclusive processes (CEP) obtained in Ref. 

using the CTEQ6L parton distribution functions. In the case of midrapidities, we also present our predictions 
considering a cutoff in the transverse momentum (pT > 6 GeV), since the CMS and ATLAS detectors have a very 
low acceptance at low pr- Our predictions for the inclusive production agree with those presented in Ref. for the 
integrated cross sections obtained without a cutoff in the transverse momentum. If the cutoff is taken into account, the 
predictions for the double J/'^ production at midrapidities are reduced by three orders of magnitude. We obtain that 
our predictions for the central diffractive production are similar to those for the central exclusive production obtained 
in Ref. [2l|. As already emphasized above, the topology of the final state of these two processes is different. While 
in central exclusive processes one has only the leading hadrons, two J/'0’s and nothing else, in the central diffractive 
case one expect to have some extra particles coming from the Pomeron remnants. Consequently, in principle, the 
experimental separation between these two processes can be performed at smaller luminosities, as those presented in 
the LHCb analysis [i^. However, it is not obvious if the double diffractive and the central exclusive mechanisms could 
be differentiated experimentally at the LHC in the next run. Unfortunately, due to the high luminosity and large 
pile-up environment the separation of the diffractive processes considering the rapidity gaps and the detection or not 
of the remnants of the Pomeron will be a hard task. The identification of diffractive processes should occur by tagging 
the intact protons in the final state using forward detectors to be installed at LHC. In this case, both contributions 
of the central diffractive and central exclusive processes for the double J/4> should be taken into account. 

Some comments about our predictions are in order. Firstly, in our estimates we are not including the feed-down from 
excited states, like J/'i’ + 'k', and relativistic corrections. Previous studies indicate that feed-down can significantly 
contribute for the double J/'^ cross section [13, while relativistic corrections decrease the magnitude of the cross 
section [Hill- Secondly, as the cross section is proportional to af and our calculations have been performed at 
leading order, our final results strongly depend on the choice of the hard scale. Clearly, the treatment of the central 
diffractive process considering the partonic cross sections at next-to-leading order is an important task for the future. 
However, the calculation of these corrections is still in progress [T5l - [l7l| . 

Finally, lets summarize our main results and conclusions. Studies of the double J/T production are expected to 
provide important insights for improving the theoretical description of the quarkonium production mechanism. In the 
last years, several groups have discussed in detail the production of this final state in inclusive processes. However, as 
recently demonstrated by the LHCb Collaboration, the study of the double J/T' production in diffractive processes is 
feasible. It has motivated the analysis of this process in the framework of the Resolved Pomeron model, which is one 
the the current models for the Pomeron. In this paper we have estimated the momentum and rapidity distributions for 
the double J/'k production in central diffractive processes at LHC energies. Our results indicate that the contribution 
of the central diffractive processes is not negligible and that its study can be useful to test the Resolved Pomeron 
model. In the future, we plan to extend our analysis for other heavy vector mesons in the final state as well as for 
single diffractive processes. 
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